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ABSTRACT 

The effect of color superconductivity on the coohng of quark stars and neutron stars with 
large quark cores is investigated. Various known and new quark-neutrino processes are studied. 
As a result, stars being in the color flavor locked (CFL) color superconducting phase cool down 
extremely fast. Quark stars with no crust cool down too rapidly in disagreement with X-ray 
data. The cooling of stars being in the Nf — 2 color superconducting (2SC) phase with a crust 
is compatible with existing X-ray data. Also the cooling history of stars with hypothetic pion 
condensate nuclei and a crust does not contradict the data. 

Subject headings: stars: evolution, stars: neutron, pulsars: general 



1. Introduction 

The interiors of compact stars have been dis- 
cussed as systems where high-density phases of 
strongly interacting matter do occur in nature, see 
(Glendenning 1996) and (Weber 1999). The conse- 
quences of different phase transition scenarios for 
the cooling behaviour of compact stars have been 
reviewed recently (Page 1992; Schaab et al. 1996, 
1997) in comparison with existing X-ray data. 

A particular discussion has been devoted to 
the idea that a strong attraction in three flavor 
wds-quark matter may allow for the existence of 
super-dense anomalous nuclei and strange quark 
stars (Bodmer 1971; Witten 1984; de Rujula and 
Glashow 1984). Thereby, in dependence on the 
value of the bag constant B different possible types 
of stars were discussed: ordinary neutron stars 
(NS) without any quark core, neutron stars with 
quark matter present only in their deep interiors 
(for somewhat intermediate values of B), neutron 
stars with a large quark core (QCNS) and a crust 
typical for neutron stars, and quark stars (QS) 
with a tiny crust of normal matter and with no 
crust (both for low B values). By QCNS we un- 



derstand compact stars in which the hadronic shell 
is rather narrow in the sense that it does not es- 
sentially affect the cooling which is mainly due to 
the quark core. However, this hadronic shell plays 
the role of an insulating layer between quark mat- 
ter and the normal crust. By QS we mean com- 
pact stars for which the hadronic shell is absent 
what allows only for tiny crusts with maximum 
densities below the neutron drip density (~ 10^^ 
g/cm^) and masses Mcr J; 10~^Mq (Alcock et al. 
1986; Horvath et al. 1991). Therefore, we suppose 
that the difference between the models of QS and 
QCNS regarding their cooling behaviour is only 
in the thickness of the crust. This results in com- 
pletely different relations between internal (T) and 
surface (Tg) temperatures for QS and QCNS. 

Recent works (Alford et al. 1998; Rapp et al. 
1998; Schafer 1998; Carter and Diakonov 1999; 
Rapp et al. 1999; Blaschke and Roberts 1998; 
Bloch et al. 1999; Pisarski and Rischke 1999; Al- 
ford et al. 1999; Schafer and Wilczck 1999a; Alford 
et al. 1999; Schafer and Wilczek 1999b) demon- 
strate the possibility of diquark condensates char- 
acterized by large pairing gaps (A ~ 100 MeV) 
in quark cores of some neutron stars and in QS 
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and discuss different possible phases of quark mat- 
ter. Large gaps arise for quark-quark interactions 
motivated by instantons (Diakonov et al. 1996; 
Carter and Diakonov 1999; Rapp et al. 1999) and 
by nonperturbative gluon propagators (Blaschke 
and Roberts 1998; Bloch et al. 1999; Pisarski and 
Rischkc 1999). 

To be specific in our predictions we will con- 
sider models of the canonical QS and QCNS of 1.4 
solar masses (1.4 M©) at a constant density. The 
constant density profile is actually a very good ap- 
proximation for QS of the mass M < IAMq, see 
Alcock et al. (1986). We consider the model of 
QCNS with a crust of the mass ~ 10" ^Mg, 
the model of QS with a tiny crust mass (Mcr ^ 
and the model of QS with no crust. For 
QCNS we shall use the same T^/T ratio as for 
ordinary NS, see (Tsuruta 1979; Maxwell 1979; 
Shapiro and Teukolsky 1983), whereas for QS we 
use somewhat larger values for this ratio, namely 
T, = 5 • 10~^T for a tiny crust (Horvath et al. 
1991) and T, = T for a negligible crust (Piz- 
zochero 1991). In the latter case, however, we as- 
sume the existence of black body photon radiation 
from the surface as for the cases of more extended 
crusts. We will estimate the cooling of QS and 
QCNS first in absence of color superconductivity 
and then in presence of color superconductivity 
for small quark gaps (A ~ 0.1... 1 MeV), as sug- 
gested by Bailin and Love (1984) and for large 
gaps (A ^ 100 MeV) as obtained in Refs. (Alford 
et al. 1998; Rapp et al. 1998; Schafer 1998; Carter 
and Diakonov 1999; Rapp et al. 1999; Blaschke 
and Roberts 1998; Pisarski and Rischke 1999; Al- 
ford et al. 1999; Schafer and Wilczek 1999a; Alford 
et al. 1999; Schafer and Wilczek 1999b). In the 
latter case wc will consider two phases: the color- 
flavor- locked wds-phase (Alford et al. 1999; Schafer 
and Wilczek 1999a; Alford et al. 1999; Schafer and 
Wilczek 1999b) and the Nf = 2 color supercon- 
ducting phase (Alford et al. 1998; Rapp et al. 1998; 
Schafer 1998; Blaschke and Roberts 1998; Carter 
and Diakonov 1999; Rapp et al. 1999; Bloch et al. 
1999; Pisarski and Rischke 1999), in which the 
s-quark is absent and the urf-diquark condensate 
selects a direction in color space whereby the color 
charge has to be compensated by the remaining 
unpaired quarks. 

Finally, we want to discuss the question whether 
the hypothesis of a color superconducting quark 



matter phase in compact star interiors is compat- 
ible with existing X-ray data. 

2. Normal quark matter 

A detailed discussion of the neutrino emissivity 
of quark matter has first been given by Iwamoto 
(1982) where the possibility of color superconduc- 
tivity has not been discussed. In this work the 
quark direct Urea reactions (QDU) d ueD and 
ue du were suggested as the most efficient pro- 
cesses. Their emissivities were estimated as ^ 

^QDU ^ 1026^^ (^^^ yXlZ ^^-3 g^^-l^ 

where at baryon densities pb — 2po the strong 
coupling constant is Qg « 1 and decreases loga- 
rithmically at still higher densities (Kisslinger and 
Morley 1976). The nuclear saturation density is 
po = 0.17fm~^, Ye = pe/pb is the electron fraction, 
and Tg is the temperature in units of 10^ K. The 
larger the density of the uds-system, the smaller 
is its electron fraction. For a density pb ~ 3po one 
can expect a rather low electron fraction of strange 
quark matter Ye ^10~^ (Glendenning 1996) and 
eq. (1) yields e^^^ ~ 10^^ T| erg cm~^ sec~^, see 
(Duncan 1983; Horvath et al. 1991). We did not 
yet discuss the strange quark contribution given by 
the direct Urea processes s ueu and ue sv. 
Although these processes can occur, their contri- 
bution is suppressed compared to the correspond- 
ing wd-reactions (Duncan 1983) by an extra factor 
sin^^c ^ 10~^, where Oc is the Cabibbo angle. 

If for somewhat larger density the electron frac- 
tion was too small {Ye < Yec — VSTrnie/ {S<Xs^^ pb) < 

2 • 10~*, for as — 0.7 and pb — 5po, rUe is 
the electron mass), then all the QDU processes 
would be completely switched off (Duncan 1983) 
and the neutrino emission would be governed by 
two-quark reactions like the quark modified Urea 
(QMU) dq uqev and the quark bremsstrahlung 
(QB) processes qiq2 qiq2i'i'- The emissivities 
of the QMU and QB processes were estimated as 
(Iwamoto 1982) 

gQMU ^ eQB(g(.r-) ^ 1020 ySg^g ^^-3 (2) 

^We note that the numerical factors in the estimates below 
axe valid only within an order of magnitude. This is suffi- 
cient for the qualitative compajrison of cooling scenarios we 
present in this work. 



2 



The latter estimate of QB cmissivity is done in the 
suggestion that the exchanged gluon is screened. If 
one suggests that quarks are coupled by transverse 
non-scrcencd gluons then one gets (Price 1980) 
eQ^(unscr) ~ lO^^ r| erg cm'^ sec'^ . With 
a nonperturbative gluon exchange (Blaschke and 
Roberts 1998; Bloch ct al. 1999) we expect that 
the estimate e^^ (scr) is more appropriate than the 
one given by e^^(unscr). Therefore, in evaluat- 
ing the cmissivity of QB processes wc will use (2). 
Other neutrino processes (like the plasmon decay 
7pi — > ee~^ vv which goes via coupling to inter- 
mediate electron-electron hole states, and the cor- 
responding color plasmon decay ffpi ^ ^ vv 
which goes via coupling of the gluon to quark- 
quark hole states (Iwamoto 1982), sec Fig. 4) have 
much smaller emissivities in the normal quark 
matter phase under consideration and can be ne- 
glected. 

Among the processes in the crust, the electron 
bremsstrahlung on nuclei gives the largest contri- 
bution to the cmissivity as estimated in (Iwamoto 
1982) 

C ~ 10^^ n erg cm-3 sec"! . (3) 

This contribution can be neglected for QS due to 
a tiny mass of the QS crust M^r ^ \^~^Mq. 

Besides, one should add the photon contribu- 
tion 

^ 2 • 10i« (y^)' erg cm'^ sec'S (4) 

where T57 is the surface temperature in units of 
10'^ K, see (Shapiro and Teukolsky 1983). This 
process becomes the dominant one for QS at es- 
sentially shorter times than for the QCNS due to 
the higher TgjT ratios for the former. 

Internal and surface temperatures are related 
by a coefficient determined by the scattering pro- 
cesses occurring in the outer region, where the 
electrons become non-degenerate. For NS with 
rather thick crust, an appropriate fit to numerical 
calculations (Tsuruta 1979) is given by a simple 
formula (Shapiro and Teukolsky 1983) 

= (10 T)2/3 , (5) 

where Tg and T both are measured in units of K. 
We shall use this expression when dealing with 



QCNS. A rough estimate of (5) yields Tg = a x 
10~^T with a ~ 0.2 — 2 in dependence on the value 
of the internal temperature varying in the interval 
10^" . . . 10^ K of our interest. In the QS the crust 
is much more thin than in the NS and the ratio 
T5/T should be significantly larger. Therefore, de- 
pending on the thickness of the crust we shall use 
two estimates for QS scenarios: Tg = 5 x 10~^T 
for a tiny crust (Horvath et al. 1991), and Tg ~ T 
for negligible crust (Pizzochero 1991). 

In order to compute the cooling history of the 
star we still need the specific heat of the electron, 
photon, gluon and quark sub-systems. In accor- 
dance with the estimates of Iwamoto (1982), Hor- 
vath et al. (1991) we have 

ce ^ 0.6 • 10^0 (^^^ Tg erg cm'^ K-\ (6) 

c-y ~ 0.3 • 10^^^ r| erg cm^^ , (7) 
Cg ~ 0.3 • 10" NgTi erg cm'^ K'^ , (8) 

c, ^ 10^1 (^^) Tg erg cm'^ R-^ , (9) 

where Ng is the number of different color states 
of massless gluons. The very small contribiition 
to the specific heat of the crust can be neglected 
(Lattimer et al. 1994). The cooling equation reads 

i=g,e,7,3 j=QDU,QMU,QB,cr 

where the summation is over all contributions to 
the specific heats and emissivities as discussed 
above. The evolution at large times on which we 
are focusing our interest here is insensitive to the 
assumed value of the initial temperature Tq. We 
checked it rising different values of initial tempera- 
ture. To be specific we choose the value To^g = 10 
as a typical initial temperature for proto-neutron 
stars. 

In Figs. 1-3 we show the cooling history of 
QCNS with standard thickness of the NS crust 
(when internal and surface temperatures are re- 
lated by Eq. (5)), QS with a tiny crust (Tg = 
5 • lO-^T) and QS with negligible crust (T, = T), 
respectively. Solid curves are for the matter sug- 
gested to be in the normal state, i. e. in the ab- 
sence of color superconductivity. Different groups 
of data points are taken from Table 3 of Ref. 
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(Schaab et al. 1999) where the notations are ex- 
plained. In the lower panels of Figs. 1-3, we show 
the results for Yg > Y^c taking Ye = 10"^, as = 1, 
p = 3pQ. This is a representative set of parame- 
ters for which the QDU processes contribute to the 
cooling. We see that two low-temperature pulsars 
can be explained as QCNS being in normal state 
(Fig. 1, thick solid curve, lower panel) and many 
observations can be interpreted as QS in normal 
state with a tiny crust (Fig. 2, lower panel). 

The upper panels of Figs. 1-3 demonstrate the 
cooling history of QCNS and QS for Ye < Fee, 
namely for Ye = 0, a., = 0.7, and p = 5po- QCNS 
being in normal state (Fig. 1, thick solid curve, 
upper panel) cool down rather slowly but are still 
in agreement with the data for a few pulsars. For 
QS with a tiny crust (Fig. 2, Te = 5 x IQ-^T) we 
get a nice fit of many data points. 

In both cases Yg < Ygc and Yg > Ygc, see Fig. 
3, QS with negligible crust cool down too fast in 
disagreement with the X-ray data. 

3. Color superconductivity 

In the standard scenario of NS cooling the in- 
clusion of nucleon pairing suppresses the emissiv- 
ity resulting in a more moderate cooling. Now, 
considering QS and QCNS we will show that we 
deal with the opposite case. 

Due to the pairing, the emissivity of QDU pro- 
cesses is suppressed by a factor exp(— A/T) and 
the emissivities of QMU and QB processes are 
suppressed by a factor exp(— 2A/T) for T < Tc- 
Thereby in our calculations we will use expression 
(1) for QDU suppressing the rate by exp(— A/T) 
and expressions (2) for QMU and QB suppress- 
ing the rates by exp(-2A/T) for T < Te. We 
also observe that plasmon and color plasmon de- 
cay processes are switched off in the superconduct- 
ing phase when the photons and the gluons ac- 
quire masses due to the Higgs effect, as it happens 
for photons in usual superconductors. Voskre- 
sensky et al. (1998) however demonstrated that 
in superconducting matter there appears a new 
neutrino neutral current process analogous to the 
plasmon decay but now due to a massive photon 
decay. Its emissivity is suppressed by the factor 
exp(— m-y/T) rather than by exp(— A/T), as for 
direct Urea processes, or by exp(— 2A/T), as for 
modified Urea and corresponding bremsstrahlung 



processes. This results in a large contribution for 
small but finite values of m-y. Naively, one could 
expect that in a color superconductor the squared 
photon mass is proportional to the fine structure 
constant a = 1/137 as it is in ordinary supercon- 
ductors. Then it would be much smaller than the 
squared gluon mass since the latter quantity in 
the color superconducting phase has to be propor- 
tional to the corresponding strong coupling con- 
stant as- In reality, due to the common gauge 
transformation for electromagnetic and color fields 
one deals with mixed photon-gluon excitations. 
We demonstrate this using the expression for the 
free energy density of the color superconducting 
phase (Bailin and Love 1984), 

/ = fn + ad*d+^b{d*d)^ 

+ c{V + iqA+'-^)d*{V-iqA^'-^)d 

{rotlf , (rotj3)^ 
^ 87r ^ 87r ' ^ ' 

where /„ is the normal part of the free energy 
density, a = jjpFqt/n^, t = {T - Te)/Te < 0, 

b = 7C(3)MPFg/(87r4T2), C = plgb/{&p?), II ~ PFq 

is chemical potential of ultra-relativistic quark, 
q = is the electric charge of a wd-pair, 

a = 1/137. We have introduced an interaction 
with two gauge fields and _B^. is the elec- 
tromagnetic field and relates to the gluons. For 
simplicity we consider only fluctuations of space- 
like components of the fields and assume the B^^ 
field to be an Abelian field. Variation of (11) 
with respect to the fields gives the correspond- 
ing equations of motion. Taking d = da + d' , 
where do = \J —ajb is the order parameter and 
d' is the fluctuating field, we linearize the equa- 
tions of motion for the fluctuating fields d',A,B. 
Solving these equations in the Fourier representa- 
tion, we get three branches of the spectrum. The 
branch lu''^ = + 2\a\ corresponds to fluctuations 
of the order parameter characterized by a large 
mass rud = -\/^2t/xpp^/7r m„ = 140 MeV. The 
branch 

uj^ = q^+mlg (12) 

describes a massive photon-gluon excitation with a 
mass rn^ g = 87rc(a -t- 3as)do/9. The extra branch 
= g ^ describes a massless mixed photon-gluon 
Goldstone excitation. 
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Thus in difference with a usual proton super- 
conducting phase of NS where the photon has 
rather low mass m-y = doVSnca ~ 4 MeV for 
fi ~ 400 McV, in the color superconductor wc, 
probably, deal with a much more massive mixed 
photon-gluon excitation (with a being replaced by 
a + 3a s) and with the corresponding Goldstone 
boson^. 

The Goldstone boson does not contribute to the 
mentioned photon-gluon decay process, whereas 
the massive excitation does. Now, armed with 
an expression for the photon-gluon mass we may 
estimate the emissivity of the corresponding pro- 
cesses (7 — 5) ee~^ + qq~^ — > vP, where 
and are the electron hole and the quark hole, 
respectively, sec Fig. 4. Using the result for 
7 ee~^ +pp~^ vv (Voskresensky et al. 1998) 
we easily get 




X e:iq){—m^^g/T) erg cm ^ sec ^,(13) 



for T < m-y^g, and using the condition A ^ /u. As 
we see, the emissivity of this process is strongly 
suppressed for the values rrij^g ~ 70 MeV following 
from eq. (12). Also the specific heat of this mixed 
photon-gluon excitation is suppressed by the same 
exponential factor exp(— m-y,g/T). For the Gold- 
stone excitation the contribution to the specific 
heat is given by Eq. (7). 

For the quark specific heat at T < Tc wc use 
an expression similar to the one which applies for 
the case of nuclcon pairing (Miihlschlegel 1959; 
Maxwell 1979; Horvath et al. 1991), i.e. 

Csg = 3.2 Cg in/T) exp(-A/T) 

X [2.5 - 1.7 T/n + 3.6 {T/Tcf] , (14) 

where Tc is related to A as A = 1.76 Tc for the case 
of small gaps. For CFL and 2SC phases we will use 
Tc ci: 0.4 A . Actually, in the latter case the rela- 
tion between Tc and A is unsettled. However, one 
believes that the cocfBcicnt in standard BCS for- 
mula Tc ~ 0.57 A is appreciably reduced as impact 

^However, the penetration depth of the external magnetie 
field is associated with the above mentioned small photon 
mass rather than with the massive photon-gluon excitation 
or with the massless Goldstone excitation, see (Blaschke et 
al. 1999). 



of instanton-anti-instanton molecules (Rapp et al. 
1999). The mentioned uncertainty in the value of 
Tc for CFL and 2SC phases does not significantly 
affect the cooling curves since the dominant effect 
comes from the exponential factor where A enters 
rather than Tc- 

Now, armed with all necessary expressions we 
may estimate the cooling of QS and QCNS being 
in uds- or 2SC phases (except for the crust). 

3.1. Cooling of different wds-phases: uds- 
smedl gaps, CFL (Fg > Yec), and 
CFL {Yc < Ycc) 

Wc select the following values of the pairing 
gaps: small gaps A = 0.1... 1 MeV as sug- 
gested to occur in the color superconducting re- 
gion by Bailin and Love (1984) and a large gap 
A ~ 50 MeV, as suggested for CFL phase in 
recent works (Alford et al. 1998; Rapp et al. 
1999). Horvath et al. (1991) have discussed the 
cooling of superconducting QS and QCNS for 
very small gaps with critical temperatures Teg = 
0.1,0.5,1. However, for QDU processes the sup- 
pression factor exp(— 2A/T) has been used rather 
than exp(— A/T). The value of the critical tem- 
perature TcQ = 0.1 seems to be quite small, there- 
fore we use A = 0.1 ... 1 MeV for the case of small 
gaps. 

We calculate the cooling history of QS and 
QCNS using eq. (10), where now summation over 
j implies summation of emissivities of QDU, eq. 
(1), suppressed by exp(-A/r), QMU and QB, 
eq.(2), suppressed by exp(— 2A/T), emissivity of 
the crust, eq. (3), and emissivity of photon-gluon 
decay, eq. (13). Summation over ? implies summa- 
tion of the quark contribution evaluated according 
to eq. (14), the electron contribution, eq. (6), 
the massless photon-ghion Goldstone contribution 
which coincides with that given by eq. (7) and the 
contributions of massive gluons given by eq. (8) 
suppressed by exp(— m^_g/r) and thus being very 
tiny. The contribution of the crust to the specific 
heat is negligible. Also in the CFL phase exist 9 
hadronic quasi- Goldstone modes. Although their 
masses are not known we may roughly estimate 
them as rrih > rUg, where rriq is the bare quark 
mass with a minimum value of ~ 5 MeV. With 
these masses the contribution of hadronic quasi- 
Goldstone modes is also very small at tempera- 
tures of our interest and can be neglected. 
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The dotted curves in Figs. 1-3 demonstrate 
the coohng history of QCNS and QS for the gap 
A = 0.1 MeV, whereas the dashed lines corre- 
spond to the cooling of the CFL phase for A = 50 
MeV. All thin dotted and dashed lines corre- 
spond to the case when the process of massive 
mixed photon-gluon decay is artificially excluded 
whereas the corresponding thick lines represent 
the cooling history when this process is taken into 
account according to Eq. (13). This new pro- 
cess essentially influences on the early stage of the 
cooling although the mass of the mixed photon- 
gluon excitation was supposed to be very high 
(m^,g = 70 MeV for Yg = 10~^ > l^ec, lower panel, 
and m-y^g = 60 MeV for Ye = 0, upper panel). This 
is due to to a big numerical factor in eq. (13). 

In all the cases we obtain very rapid cooling 
in disagreement with the data. Particularly rapid 
cooling occurs for the CFL phase. In the latter 
case contributions of the QDU, QMU and QB pro- 
cesses to the emissivity arc suppressed as well as 
the quark contribution to the specific heat. The 
rate is governed by the photon emissivity from the 
surface. For Ye > Yec, the specific heat is deter- 
mined by the electrons. As the consequence of this 
reduction of the specific heat we get a very rapid 
cooling of the CFL (Ye > Yec) phase, see lower 
panel of Figs. 1-3. For Ye = (upper panel of 
Figs. 1-3) there are no electrons and the specific 
heat is determined by a very small contribution 
of the Goldstone mode given by eq. (7), so that 
both QCNS and QS cool down even faster than 
for the case Ye = 10~^ > Fee- In both Ye > Yec 
and Ye < Yec cases the cooling time of the CFL 
phase is extremely small. This means that in re- 
ality the cooling is governed by the heat transport 
in the thin crust (Pizzochero 1991), which we did 
not take into account. 

Thus we sec that QCNS and QS, if present 
among the objects measured in X-rays, can't be 
in the CFL phase. The cooling of this phase is 
so rapid that one might expect problems not only 
for the models of QS and QCNS but also for the 
models of NS with quark cores consisting of the 
CFL phase only in deep interiors. 

3.2. Cooling of the 2SC phase 

This phase is probably more reliable for QCNS 
rather than for QS since the CFL phase is ener- 
getically favorable in the latter case. The 2SC 



phase is characterized by large gaps, A ^ 100 MeV 
(Alford et al. 1998; Rapp et al. 1998, 1999). To 
be specific we suppose that blue-green and green- 
blue ud-quarks are paired, whereas red u- and d- 
quarks {ur, dr) remain unpaired. This has the con- 
sequence that the QDU processes on the red (un- 
paired) quarks, as dr UrcD, as well as QMU, 
drQr — * Urqre.v, and QB, qirCtir QirQ2ri^t/, are 
not blocked by the gaps whereas other processes 
involving paired quarks arc blocked out by large 
diquark gaps. The QDU process on red quarks 
occurs in the Yg > Yec case only. Its emissivity is 
given by 

erg cm~^ sec~^. (15) 

The extra suppression factor of the rate (1) comes 
from the fact that the number of available un- 
paired color states is reduced. 

QMU and QB processes on red quarks are also 
rather efficient. Although there is no one-gluon 
exchange between dr — dr, the QMU and QB pro- 
cesses may go via a residual quark-quark inter- 
action, e.g, via two-gluon exchange. We roughly 
estimate the corresponding emissivities as 

~ lO^^Tl ergcm"^ sec'^ (16) 

In the 2SC (Ye > Yec) phase the QDU process 
on red quarks is the dominant process and QMU 
and QB processes on red quarks are subdomi- 
nant processes whereas in the 2SC (Ye < Yec) 
phase QDU processes do not occur and QMU and 
QB processes on red quarks become the dominant 
processes. Other processes like QDU, QMU and 
QB with participation of other color and fiavor 
quarks are continued to be appreciably suppressed 
by large gaps. 

The specific heat is also changed in the 2SC 
phase since the dr and Ur contributions are not 
suppressed by a factor exp(— A/T) whereas color- 
paired ud-contributions remain to be suppressed. 
With these findings we calculate the cooling his- 
tory of QCNS and QS. The results are presented 
in Fig. 5 for YJ. = 10^^, p = 3po (thick lines), 
and Ye = 0, p = 5po (thin lines). We see that in 
both cases the cooling history of QCNS and also 
of QS with a tiny crust (T^ = 5 • lO'^T) nicely 
agrees with the X-ray data. The cooling of QS 
with negligible crust does not agree with the data. 
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4. Conclusions 

We have estimated the contributions of vari- 
ous quark processes to the emissivity. Among 
them, the new decay process of the massive mixed 
photon-gluon excitation is operating at the early 
stage of the coohng and QDU, QMU and QB pro- 
cesses on red quarks determine tlie cooling of the 
2SC phase. We discussed the cooling history of 
QS and QCNS taking into account different pos- 
sibilities: Ye > Yec and Yg < Yec, the normal 
quark phase, and various color superconducting 
phases as the "uds-phase - small gaps" suggested 
by Bailin and Love (1984), the CFL phase, and the 
2SC phase, as suggested in recent works (Alford ct 
al. 1998; Rapp et al. 1998; Schafer 1998; Alford et 
al. 1999; Rapp et al. 1999). In all the cases we see 
that QS and QCNS being in the CFL phase cool 
down extremely fast in disagreement with known 
X-ray data. Also the cooling curves for the case of 
small gaps (A = 0.1 ... 1 MeV) disagree with the 
data. 

Even if the CFL phase would be realised only 
in the deep interior region of a NS it would be 
problematic to satisfy X-ray data. In this case the 
star would radiate mostly not from the surface but 
from the CFL region due to its extremely small 
specific heat related to the Goldstone excitation. 
Thus the cooling time would be determined by the 
heat transport from exterior regions to the center 
rather than by the cooling of the hadronic shell. 

The cooling history of the QS and QCNS with 
a crust being in the normal state agrees with the 
data. 

In this respect the following remark is in or- 
der. It is now believed that quark matter below 
Tc ^ 50 McV is in the color superconducting state 
characterized by a diquark condensate with large 
energy gaps (A ~ 100 MeV) rather than being in 
the normal state or the superfluid state character- 
ized by small gaps (A < 1 MeV). If so, one could 
think that our above discussion of normal quark 
matter and of the case of a small gap has just ped- 
agogic reasoning. However, this is not really so. 
Indeed, besides the idea of abnormal strange nu- 
clei and strange stars (Bodmer 1971; Witten 1984; 
de Rujula and Glashow 1984) there is the very 
similar idea of abnormal pion condensate nuclei 
and stars with pion condensate nuclei, see (Migdal 
1971; Voskresensky 1977) and the review (Migdal 



et al. 1990), chapters 15, 16. The same relates to 
the kaon condensate objects. Pion condensate sys- 
tems cool down at about the same rate as given by 
QDU processes (for ^ 10~^) ^. Besides, they 
can be in the normal state or in the superfluid 
state characterized by very small gaps A < 0.1 
MeV. The cooling history of systems being in nor- 
mal state is described by thick solid curves on the 
lower panels of Figs. 1-3. Thus we may also 
conclude that the hypothesis of pion condensate 
nuclei-stars (being in normal A = state with a 
crust) does not contradict to the X-ray observa- 
tions. Stars with pion and kaon condensate nuclei 
being in the superfluid state with gaps A > 0.1 
MeV are ruled out as objects being observed in 
X-rays. 

The cooling history of the 2SC phase of QCNS 
and QS with a tiny crust (Tg = 5 ■ 10~'^T) agrees 
with the X-ray data. The cooling history of QS 
with no crust disagrees with X-ray data. 

Three final remarks are in order: 

(i) It is conceivable that there are more complex 
collective effects which essentially affect the spe- 
cific heat and the luminosity. E.g., we calculated 
the mixed photon-gluon spectrum in a simplified 
model of two Abelian gauge fields and concluded 
that the mass of the excitation is large, whereas 
one can't exclude that in the realistic non- Abelian 
case there exists a photon-gluon excitation of a 
small mass that could lead to very efficient cooling 
via the mixed photon-gluon decay process given by 
eq. (13). The masses of hadronic quasi-Goldstone 
modes in CFL phase should be carefully studied. 

(ii) As we mentioned above in the discussion 
of the DU process, we have neglected the contri- 
bution of strange quarks (QDU-s) relative to that 
of the light quarks. The former contribution is 
^ 10~^ times smaller than the latter in the case 
of normal matter and in the case when all diquark 
gaps are identical, the discussion given for the CFL 
phase applies. However, if the pairing gaps for 
strange diquarks are smaller than those for non- 
strange diquarks the QDU-s contribution to the 
emissivity can be essentially enhanced. Up to now 

'One should bear in mind a strong suppression of emissiv- 
ities of pion condensate processes due to nucleon-nucleon 
correlation effects (Blaschke et al. 1995) and an enhance- 
ment of the specific heat (Voskresensky and Senatorov 
1984, 1986; Migdal et al. 1990) which are often ignored 
in the cooling simulations. 
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there exist only rough estimates of the values of 
the gaps and this gives rise to large uncertainties 
in final estimates of the emissivity. Above, in or- 
der to be specific, we made calculations consider- 
ing QDU on u- and d- quarks only. The inclusion 
of QDU-s in the case when the strange diquark 
gaps can be smaller than those for non-strange di- 
quarks is simulated by varying the values of the 
gaps in a wide interval. This does not change the 
qualitative picture of the QS and QCNS cooling 
we have discussed in the present work. 

(iii) We also would like to point out that, if 
the compact object formed in the explosion of SN 
1987A was a QS or a QCNS being in the CFL 
phase {Ye < Yec), it is now so cold that it is al- 
ready impossible to observe it in soft X-rays. This 
becomes particularly interesting if continued ob- 
servation of SN 1987A would find a pulsar and 
would not observe it in X-rays. 
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fruitful discussions by M. Alford and K. Rajagopal 
after reading the draft of the paper. We also thank 
B. Friman, J. E. Horvath, E. Kolomeitsev, R. Pis- 
arski, D. Rischke, A. Sedrakian, and F. Weber for 
their discussions. One of us (DNV) is grateful for 
the hospitality extended to him during the visit 
at Rostock University and acknowledges financial 
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Fig. 1. - Cooling history of QCNS M = 1.4 Mq. 
The relation between Tg and T is given by eq. (5). 
Lower panel: Ve = 10~^, p = 3po, and as = 1; 
upper panel: = 0, p = Spo, and = 0.7. 
For the thin dotted and dashed lines the process 
of massive mixed photon-gluon decay is artifically 
excluded; the corresponding thick lines take this 
process into account. Line styles are related to 
different values of the gap A, given in the legend. 
Data points are taken from analysis Schaab et al. 
(1999). 

Fig. 2. — Same as in Fig. 1 but for QS with T, = 

5 • IQ-'^T. 

Fig. 3. — Same as in Fig. 1 but for QS with T, = 
T. 

Fig. 4. — Photon- gluon decay via intermedi- 
ate quark- quark hole and electron- electron hole 
states. 

Fig. 5. — Cooling history of QCNS and QS being 
in 2SC phase, M = 1.4 Mg, A = 100 MeV. Thick 
lines: Yg = 10~^, p = 3po, and = 1; thin lines: 
Ye = 0, p = 5poi a-nd = 0.7. Solid curves 
correspond to QCNS, the relation between Tg and 
T given by eq. (5), dashed lines correspond to QS 
with a tiny crust, T, = 5 • 10~^T, and dotted ones 
to QS without crust, T, = T. 
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g. 1: Blaschke et al. 




g. 2: Blaschke et al. 
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Fig. 4: Blaschke et al. 
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